When physics helps mathematics: calculation of the sophisticated 
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There exists a remarkable connection between the quantum mechanical Landau-Zener problem 
and purely classical-mechanical problem of a ball rolling on a Cornu spiral. This correspondence 
allows us to calculate a complicated multiple integral, a kind of multi-dimensional generalization of 
Fresnel integrals. A direct method of calculation is also considered but found to be successful only 
in some low-dimensional cases. As a byproduct of this direct method, an interesting new integral 
representation for £(2) is obtained. 
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Abstract 



INTRODUCTION 

According to Vladimir Arnold [1], mathematics can be considered as some branch of 
physics. In writing this note we have no intentions to advocate such point of view. Never- 
theless, in our opinion the following calculus problem is a hard nut to crack if only mathe- 
matical considerations are being used. The problem lies in exactly calculating the multiple 
integral of the type 

OO Si S2n— 1 

I n = J ds i J ds 2 --- J ds 2n cos Oi - si) •••cos(s2 n _i-S2 n ). (1) 



- OO 



Because of the si -H- s 2 symmetry, the n = 1 case is simple 

_^ OO OO 

Ji = — J dsi J ds 2 (coss^ cos s\ + sin s\ sins^) = — , (2) 

— OO — OO 

since its calculation involves known Fresnel integrals 



ds cos s 2 = I ds sin s 2 
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— oo — OO 

However, already for n > 2 the above symmetry is lost and things quickly become messy. 
The n = 2 and n = 3 cases lay at the borderline. They can be done with some efforts even 
though the calculations become noticeably more involved. Unfortunately, they do not admit 
an apparent generalization by using the induction method. Already for n = 4 the attempt 
to use the same methods meets difficulties. Nevertheless, we found a way to calculate such 
type of integrals by invoking some physical arguments. 



"IF YOU CANNOT SOLVE A PROBLEM, THEN THERE COULD BE AN EASIER 
PROBLEM YOU CAN SOLVE" 



Trying to follow George Polya's advice, let us consider the following system of ordinary 
differential equations 
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(3) 



2 



Such system of equations emerges naturally when one tries to describe a motion of a sphere 
S 2 rolling on the flat surface R 2 without slippage. More accurately, it describes the rolling 
of a sphere of radius R along the Cornu spiral on R 2 whose curvature k = as is proportional 



to the arc-length s [3, |4|. But what is the relation of this problem to our integral ([T])'. 



As the matrices 



M{s) 



I 









sin 



cos ■ 



sin ^ - cos ^ 



J 



do not commute for different values of s, the solution of ([3]) is given by the time ordered 
exponential 

U(s, s ) = Texp ( J M(si) dsA = 1 + J M(s x ) ds Y + 



\*0 



*0 



S Si 



J d Sl J ds 2 M{s l )M{s 2 ) + J dsi J ds 2 J ds 3 M{s x ) M{s 2 ) M{s 3 ) + ■ ■ ■ (4) 



SO SO 



SO SO SO 



so that 



( x(s) \ 



y[s) 



U(s,s ) 



y(s ) 
\4*o) J 



(5) 



Let us take a closer look at the individual terms of the above infinite series (jlj). Writing the 
matrix M(sj) in the block-form 

/ 

Mi EE M( Si ) = 



X l 

-xl 



where 



Xi 



R 



- sm f 

as 2 ! ' 

COS 



it is easy to prove by induction that 



M X M 2 ■ ■ • M; 



2n+l 



{-^-TXlXl ■ ■ ■ X2n-lXl n X2n+l 



\ (-i) n+1 XlX2 ■ ■ ■ xL-lX2nXln 



+ 1 



and 



/ 



M X M 2 ■ ■ ■ M- 



2n 



{-^) n X\Xl---X2n-\Xln 



\ 



(-^) n XlX2---xIn-lX2n ) 



But 



T 1 a I 2 2 \ 

Xi Xi+i = cos -{Si - s i+1J 



Besides, then our calculations is supplemented by the following initial and boundary condi- 
tions 

so = — oo, s = oo, x(— oo) = y(—oo) = 0, z(— oo) = 1, (6) 
in view of these relations and taking into account (jSJ) and (j3j), we finally obtain 

oo ,-Qn ~ J s y 

Z(00) = 1 + Xj ffin 7 rfSl _/ rfs 2--- J ^S 2 „ /„(Si,S 2 , • • • ,S 2 „), (7) 

where 



-oo — oo 



r , . a /2 2\ a /2 2\ J> 

/ n (Si, S 2 , . . . , S 2n ) = COS - S 2 ) COS -(S 3 - S 4 ) • • • COS -(San-1 - S 2n) 



After rescaling 



Si ^ J -Si, (8) 



a 

our original integral ([1]) indeed shows up in the series (jTJ): 

oo / 9 \ n 

*«») = i + E(-i)-(^) 

Thus, if we can find -2(00), we can find I n as well! But how can we find z{oo)l 



(9) 



SPINORIZATION AND THE HOPF MAP 

Now we shall follow the advice of Jacques Hadamard: "The shortest path between two 
islands of truths in the real domain passes through the complex plane". By adopting his 
suggestion to our case it is convenient at this point to use two complex variables a and b 
instead of three real variables x, y and z through relations [5j 

x = ab* + ba*, y = i(ab* — ba*), z = aa* — bb* . (10) 

Notice: introduction of two complex variables is equivalent to looking at solution for our 
problem in C 2 ! Furthermore, in view of the initial conditions and ()3]), the variables x,y and z 
are constrained to unit sphere S 2 . This causes the variables a and b to be constrained to S 3 , 
that is to obey the equation \a\ 2 + \b\ 2 = 1. See Ref. jj| for more details. Interestingly enough, 



4 



under these conditions Eq{10] describes the Hopf map S 3 -> S 2 Q. From Ref.l5j it can be 
seen that the complex variables a and b must satisfy the following system of differential 
equations 



.d_ 

ds 





1 / 






(;) 


2R 1 





ias 2 /2 



-ias 2 /2 








(11) 



if the real variables x, y and z to satisfy ([3]). Now we can again formally solve (ITTj) by using 
the time-ordered exponential series. This time, however, the solution is known and it was 

n 

obtained by Rojo in |2J. We just shortly repeat it to ensure the continuity of our exposition. 
The solution is formally given by 




U(s, s 




where 



with 



S S Si 

U(s, s ) = 1 - i J H(si) d Sl + (-if J dsi J ds 2 H( Sl ) H(s 2 ) + 



(12) 



(13) 



so 



so 



so 



His) = - 



2R 



Q e -ias 2 /2 
3 ias 2 /2 q 



The imposed initial conditions ()6]) are now translated into 

a(— 00) = 1, b(— 00) = 0. 



(14) 



Since the unitary evolution ( TT2|) conserves the norm aa* + bb*, from (I14p we obtain back the 
equation for 3-sphere, that is aa* + bb* = 1, valid for any "time" s. This result implies that 



z(oo) = |a(oo)| 2 - |6(oo)| 2 = 2|a(oo)| 2 - 1. 



(15) 



Evidently, we need only to calculate a(oo) to obtain z(oo). 

By examining the product of H(si) matrices, we observe that for the odd number of 
multipliers the matrix product does not have non-zero diagonal terms and, hence, does not 
contribute to a(oo) thanks to the initial conditions (|14p . The remaining terms with even 
number of multipliers have easily calculable non-zero diagonal elements so that we get 

2 OO Sl 



a 00 



1 + (^) / d Sl e-^ 2 J ds 2 e^ 2 + • • • . 
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After rescaling flSJ) this expression acquires the form 



(°c) = ! + £(-!)" 



n=l 



2aR 2 



Jni 



where 



si 



J n 



ds 1 e iS i / ds 2 e™ 2 / ds 3 e 183 ■ ■ 



ds 2n e ls ^. 



(16) 



(17) 



In contrast to I n , the multiple integral J n is doable. It can be calculated as follows [2j. First, 
we write 



J n = J dsx--- J ds 2n 6{s x - s 2 ) ■ ■ ■0{s 2n -i - s 2n ) e n (sl, . . . ,s 



2 ) 

2n)i 



with 



I 2 2 \ _ J_7 2 _ 2 _l_ 2 _ 2 _l_ _l_ 2 _^ 2 H 

e n(,'5 1 , . . . , s 2n ) — exp I l[S 1 S 2 + s 3 s 4 -+- ■ ■ • -+- s 2n _ 1 — s 2n )\- 



Then, we use the integral representation for the Heaviside step function 



9(s) 



2m 



dbj 



oj — le 



1, if s > 0, 
|, if s = 0, 
0, if s < 0. 



(18) 



Using the result, we can perform the integrals over dsi. For this purpose we sequentially 
complete the squares, e.g. 

cu 2 -uji\ 2 (oj 2 - UJi) 2 



1 s 2 [u 2 - Ui) + s 2 
and then evaluate the Gaussian integrals 



S2 + 



4 



e ±is2 ds 



(cos s ± isins ) ds — \ — (1 ± i) 



As a result, we finally obtain 



7 

Jn = 7T 



duii 
2tH 



dui- 



2n~l 



2m 



cxp 



1 1 [£J 2 (Wi - W 3 ) + W 4 (^3 - W 5 ) H h W 2n _2(w 2 n-3 - ^2n-l)]} 



(wi-ze)(w 2 -ie)---(ct;2 n _i-ze) ' ^ 

As it can be seen, the terms quadratic in Ui are all canceled thanks to the alternating signs 
in exponents in ( 1T7|) . It is this feature that distinguishes, as we shall demonstrate in the next 
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section, the calculation of J n from the calculation of /„ and makes the integral J n solvable. 
For this purpose we rescale even-index variables u 2 i — > 2u 2 i in f fl9|) and perform integrals in 
these variables taking into account (TT8|) . The result is: 



oo oo 

duii f dec?, 



2iri J 2m' 

oo — oo 



duj 2n -i - W 3 )#(w 3 -u s )--- d{u 2n -3 - UJ 2 



n-lj 



or 



2m (ui - ie)(u 3 - ie) ■ ■ ■ (uj 2n -i - ie) 

oo 

aw i f du3 f du 2n -\ 



2m J 2iri 

oo — oo 



7T 



(20) 



27ri (wi - ie)(u> 3 - ie) • • ■ (w 2 n-i - ie) ' 

It is the symmetry of the integrand in (1201) that makes the calculation of fl20l as easy as the 
calculation of I\ from which our story had begun. In the present case we obtain 



7 

Jn = 7T — 



duj\ 1 



27ri wi — 



7r n 1 /n\ n 

nW®' = n(i) ■ 



Then ffT6|) shows that 



and from flT5|) we get 



a( °° ) = eXP (~i^) : 



(21) 



(22) 



z(oo) = 2 exp 



7T 



2ai? 2 



°° 2 / 7T V 

1 = 1 + £(-i>"! (d* 



n=l 



Comparing with (|9]), we finally obtain the desired expression for the integral I n : 

2 (ix^ n 



n\ V4 



(23) 



DIRECT CALCULATION OF J 2 



If you are still unhappy by our use of indirect methods of calculation of deceptively simply 
looking integral ([I]), here we discuss some features of the direct method. Unfortunately, as 
far as we can see, it works well only for small values of n. 

Indeed, let us write 

oo oo oo oo 

h= J ds\ J ds 2 J ds 3 J ds^ 04 (si, s 2 , s 3 , s 4 ) cos(sj-sl) cos (s^ - S4), (24) 
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where 



4 (si, s 2 , s 3 , s 4 ) = 0(si - s 2 ) 6>(s 2 - s 3 ) 6>(s 3 - s 4 ) 



(25) 



Then, as in the previous section, we use the integral representation for the step function and 
notice that 



OO CO 



OO CO 



J d Sl J ds 2 e iu)lSl e- iS2{u)1 -" 2) cos {si -s 2 2 ) = )- j ds x J 



-CO —CO 



-co -co 



ex v[- l T + l [ Sl + T) +l — 4~ 



M S 2 + 



U)i - u 2 



\+ 



exp I z— — % \ s\ 



Ul\ 2 _ , (^1 - ^2) 



2 J 



71 COS 



+ M s 2 - 



U)i - u 2 



4 



cu 2 (2wi - w 2 ) 



(26) 



Similar calculations are done for integrals over ds 3 and <is 4 . As a result, after rescaling 
uj 2 —7- 2w 2 , we end up with the result: 

duii °f doj 2 °f dco 3 cos [co 2 (ui - co 2 )] cos [w 2 (w 3 - co 2 )] 

[27) 



I 2 = 7T ! , , 

2ixi J 2ixi J 2ixi (tui — ie) (co 2 — ie) (w 3 — ie) 

— mo — CO — CO 

The integrals over du\ and du; 3 can be easily calculated. Indeed, 
duj\ cos (u\U)2 — w 2 ) 1 f duj\ 



2i\i oj 1 — ie 



2ni 



_T e -^2 -If- _ e w 2 

uji — ie U\ — ie 



Therefore, 



h = — 



7T 2 f dt0 2 



9(u 2 )e-^ + #(-w 2 )e^ 2 



2ni 



u> 2 — ie 



Now we use the well known result 



— - — = P— + iir6(u)) 
uj — ie lo 

to split the integral (|28p into the principal value and the 5-function parts: 

7T 2 

h = ~^(^2P + I25) ■ 
Of course, the 5-function part is obtained instantly 



(28) 



(29) 



1 



'25 
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As for the principal value part, we have 



lop = lim 



lim 

2ni e-+o 



OO o ■ 2 



W 2 



du 2 + 



Lu'2 



dujo 



duj<2 



1 /■ sin (2^2) 



duo — — — 



0O 2 J U! 2 

1 7sm[(V2uj 2 ) 2 



2tt J (^2oo 2 ) 



d[(V2u 2 f 



In making the last step we have used the Dirichlet integral 



sm uj 71 

du = — . 

u 2 



Putting all terms together, we obtain finally 



7T 

16' 



(30) 



Which is a special case of (J23J), as expected. 

Can we think about the general case (n > 2) by acting in the manner just described? By 
repeating the above steps when n > 2, we obtain 



7 du 2 




00 


J 2m 

—00 


J 2m 

—00 


■/ 

—00 



duJ 2n -2 fn[U 2 , W4, . . . , U 2n - 2 ) 



2m (u 2 - ie)(u 4 - ie) ■ ■ ■ (u 2n ^ 2 - ie) 



(31) 



where 



f n (uj 2 , LO4, ■ ■ ■ , ^2n-2) = 0(0, UJ 2 )(P(UJ 2 , W 4 ) • • • 0(w 2 n-4, ^2n-2)0(^2n-2, 0), (32) 



with 



(wi,w 2 ) = - w a ) e- 4{w i-^ 2) + #(w 2 - wi) e^-"?). 



(33) 



Evidently, for n > 2 things begin to look rather inconclusive and the above direct method 
needs some fresh input in order to be brought to completion. 



DIRECT CALCULATION OF J 3 

Now let us calculate 



7T 



2 1 '■ 



9 



where 

dx <fi(0,x)f(x) 



oo 

/ = 



and 



27rz x — ie 



f( x ) = / 7T- : ' 

2i\i y — ie 



Using the relation fl29|) in fl34"|) . we get 

~ dx <f>(0,x)f(x)-(f>(0,-x)f(-x) 



1 = 1^(0, 0)/(0) + lim / — 



.r 



However, 0(0, 0) = 1, while 

/(Q) = Jdy_ 0(O,y)0(y,O) = ±^ (0j Q) + Um J dy_ 2 (O, y) - 2 (O, -y) 



27ri y — ie 2 e->o 7 2iri y 

— oc 6 

and since 

2 (O,y)-0 2 (O,-y)=e- 2 ^-e 2 ^ 2 , 

when y > e > 0, we get 

1 1 7sin(2y 2 ) 1 



2 Ti J y 4 
o w 



and, therefore, 



/ = ^ + K, 



where 

, /" dx e~ lx fix) — e lx f(—x) 

K = hm / 

t^oJ 2ni x 



Now, using again ( 1291) . we have 



f(x) = ^(x, 0)0(0, 0) + hm " ^ 0) - 0) 



2" v ' /rv ' ' e->o7 2?n y 
and since in ( |37l) x > 0, we get 

= i e _« 2 lim 7 rfy e(x - y)e- 2 + % - x)e-^ 2 -^ - e -^ 2 -^ 2 > 

M ' 2 e-W 2m y 
Analogously, 

= 1 tf Um 7 dy e-^ 2 - 2 ) - 9{x - y)e" 2 - 6(y - x)e~^~^ 

JK 1 2 e^oj 2ni y 
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In light of (EZD, ([3HD and ([39]), 

K = K 1 + K 2 , 

where 

1 f dx e -e 1 / > sm(2a; ) , 1 
Ki = = / dx = — 

2 J 2ni x 2ir J x 8 


and 

f dx f dy 9(x — y) cos (2x 2 ) + 9(y — x) cos {2y 2 ) — cos 2(x 2 — y 2 ) 

K 2 — - 
After rescaling 



x J y 2n 2 
y 



x y 
x — > — =. -u — >■ — =, 



and using 

A ,. 1 ,7/ ,.2 \ 

0(x - y) = #(x 2 - y 2 ), if x > and y > 0, 



g?x 1 <i(x 2 ) 
x 2 x 2 ' 



as well as 

cos (x — y) = [6(x — y) + 0(y — x)] cos (x — y), 
we end up with the result 

1 [ dx J dy 



I = Ko 



f dx f dy 

J — J — [cos (x — y) — cos x 



4ir 2 J x J y 




To calculate the integral 



00 . X 



~ f dx f dy 

1 = J — J — [cos (x — y) — cos x 



x J y 




we introduce a parametric integral related to it: 



00 , X 



1(a) — ( — f — [cos (x — ay) — cos x 
J x J y 
y 

Note that 1(0) = and 1(1) = J, so that 



I - } d -Mda. 

da 



However, 



— 1^1 = f — f dy s in(x — ay) = — lim / — [cos (1 — a)x — cos x 1 , 
da J x J yi a^oi x L J 
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which is the same as 
dl{a) 1 



lim 

da a e->o 



_e(l-a) 



dx f dx 

— cos x — — cos x 
x J x 



lim 



Ci{e) - Ci{e{\ - a)) 



a 



where 



Ci{x) 



dx 



cosx 



x 



stands for the integral cosine function. Using the well known series representation for this 
function 

oo (_ T 2\k 



7 ~ 0.5772 being the Euler constant, we get 

df{a) 1 



da a £->o 



lim[ln e — In (1 — a)e] 



ln(l -a) 



a 



and, therefore, 



ln(l -a) 



7T 



da = C(2) 
a 6 



We have just proved an interesting identity which seems to be a new integral representation 
for C(2): 

dx f dy 
V 

o o y 

Collecting all pieces together, we get finally / = 1/24 and 



C(2) = J — J — [cos (x — y) — cosx]. 



(41) 



vr 3 1 2 /7T\ 3 



24 3! V47 ' 



in agreement with 



CONCLUDING REMARKS 



This problem had originally aroused in the context of a remarkable correspondence be- 
tween the quantum mechanical Landau-Zener problem (known in the context of molecular 
scattering) and purely classical problem of a ball (that is 2-sphere) rolling on a Cornu spiral 
(that is on the curve in R 2 , known as Cornu spiral) recently established by Bloch and Rojo 
in 0, Q. In fact, the main ingredient of this connection - the app.ication of the Hopf ,nap 
- goes back to Feynman, Vernon and Hellwarth |5( who showed that the quantum evolution 
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of any two-level system is determined by the classical evolution (precession) of the magnetic 
dipole moment of unit strength in an effective external magnetic field. 

It is remarkable that physics helps to calculate a complicated integral ([T]). However, we 
suspect that there should be a direct method of calculation. For low values of n, we have 
provided some examples of the direct method. It is after the readers to tackle the case of 
general n. 
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